forvalues j=1/15{
	use sim`j'.dta, clear
	gen t1 = gap/gap_se
	gen t2 = lm/lm_se
	gen t3 = me/me_se
	gen r1 = abs(t1)>invnormal(0.975)
	gen r2 = abs(t2)>invnormal(0.975)
	gen r3 = abs(t3)>invnormal(0.975)	
	disp "`j'"
	tabstat _all, col(stat) stat(mean sd p1 p5 p10 q p90 p95 p99) format(%8.4f)
}
